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Wavefunctions of complex lanthanide atoms Ndl, Pml and SmI, obtained via multi-configuration 
Dirac-Fock method, are analyzed for density of states in terms of partial densities, strength functions 
{Fk{E)), number of principal components (^2(-E)) and occupancies {{ria)^) of single particle orbits 
using embedded Gaussian orthogonal ensemble of one plus two-body random matrix ensembles 
[EG0E(H-2)]. It is seen that density of states are in general multi-modal, Fk{E)'s exhibit variations 
as function of the basis states energy and ^2{E)'s show structures arising from localized states. The 
sources of these departures from EG0E(l-|-2) are investigated by examining the partial densities, 
correlations between Fk{E), (,2{E) and (ua)^ and also by studying the structure of the Hamiltonian 
matrices. These studies point out the operation of EGOE(l-|-2) but at the same time suggest that 
weak admixing between well separated configurations should be incorporated into EGOE(l-(-2) for 
more quantitative description of chaos and localization in Ndl, Pml and SmI. 

PACS numbers: 05.30.-d,05.45.Mt,32.10.-f,32.30.-r 

I. INTRODUCTION 

Lanthanide atoms exhibit complicated configuration mixing and have complex spectra [1] as consequence of par- 
tially filled high angular momentum 4/ and 5c? valence shells. The general form of the ground state configuration is 
[Xe]6s^4/™5c?", where m and n are the shell occupancies. Across the period to increases from 1 to 14, whereas n is 1 
for Ce I and Gd I which are located at the beginning and middle of the period and n is for the remaining lanthanide 
atoms. Along the period the number of configurations and mixing increases as 4/ shell occupancy increases, and 
reaches a maximum around the middle. Lanthanide atoms are appropriate to study complexities emerging from two- 
body interactions and a comparative study along the period could provide insights to the implications of increasing 
particle number to the nature of the complexity. 

Detailed studies of Ce I wave-function structure, first carried out by Flambaum and collaborators [2, 3] showed 
the smoothed (with energy) strength functions (Fk{E)) and number of principle components (NPC or ^2iE)), which 
measure chaoticity, can be understood in terms of banded random matrices (BRM) with the local level and strength 
fluctuations following the predictions of Wigner's Gaussian orthogonal ensemble (GOE) of random matrices. Later 
studies of Cummings and collaborators on Pr I (and also Ce I) [4], the element next to Ce I in the lanthanide series, 
showed results in agreement with Flambaum et al studies of Cel. Flambaum and collaborators further developed, 
using the BW form for Fk{E) and NPC in terms of the BW spreading widths which are a result of many body 
chaos, a statistical theory for the distribution of occupation numbers [5, 6], and also for electromagnetic transition 
strengths in atoms with complex (or chaotic) spectra [7]. Applying this statistical spectroscopy for complex states, 
energy-averaged cross sections are obtained in terms of sums over single electron orbitals and using this Flambaum 
et al [8] and Gribakin et al [9] explained the observed, and quite puzzling, low-energy electron recombination rate of 

Following the successful applications in nuclei [10-13], mesoscopic physics [14], in the context of quantum com- 
puters [15] and as suggested more recently for example in the dynamics of cold atoms [16], a better random matrix 
hypothesis for atoms is to consider the interaction (after subtracting the single particle field generated by the Coulomb 
Hamiltonian see Eq. (5) ahead for the Hamiltonian) to be random. This gives rise to embedded GOE of one plus 
two-body interactions [EG0E(l-|-2)]. See Sect. HI ahead for details of EG0E(l-|-2) and here it suffices to say that 
the ensemble is {iJ} ~ h{l) + A{y(2)} where h{l) is the mean-field one-body part and {T^(2)} is a GOE in 2-particle 
space with A the interaction strength. An extreme limit, with strong interaction (A oo), of EG0E(l+2) is two body 
random matrix ensemble (TBRE). Recently we demonstrated that Sm I [17], which has three more active electrons 
compared to Pr I, carries TBRE signatures which include Gaussian state densities, strength functions having Gaussian 
form, NPC as a function of energy is of Gaussian form, local level fluctuations follow GOE etc. Going beyond this 
limiting situation, EG0E(l+2) with increasing interaction strength (equivalently with increasing particle number if 
the interaction is fixed), strength functions exhibit Breit-Wigner to Gaussian transition. It is important to stress that 
EG0E(l-|-2) at the weak limit (A small) shares features with BRM, for example the BW form for strength functions 
(this explains partly the applicability of BRM for Cel and PrI). In a recent paper [18], an interpolating function, anal- 
ogous to the Brody distribution for nearest neighbor spacing distribution, representing the BW to Gaussian transition 
of the strength function is reported (see Eq. (23) ahead) and it is used to understand this transition in rare earth 
atoms (at the spectrum's center) as we go from Cel to SmI with valence electrons changing from 4 to 8. In addition, 
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it is noticed in Ndl, Pml and SmI studies that the eigcnfunctions exhibit locaHzation properties. The presence of 
locaUzed states is also evident in the study of Ce I and Pr I by Cummings and collaborators [4]. 

As shown in [19], the number of allowed 4/™6s^ — 4/™6s6p and 4/'" — 4/"'^^5(i dipole transitions of the rare 
earth atoms is large ( for example, in Pr I the number of allowed 4/^6s^ — 4/"^6s6y» dipole transitions is 7402) and it 
is difficult to study the transitions individually. However, statistical spectroscopy, based on TBRE and many-body 
chaos, provides a means to understand the smoothed part of these transitions (with fluctuations following GOE). 
In this paper we consider Ndl, Pml and SmI with 6,7 and 8 valence electrons and examine in some detail: (i) the 
extent to which EG0E(H-2) applies for these atoms; (ii) localization properties of wavefunctions; (iii) the structure 
of the Hamiltonian matrix. For (i) and (ii), wavefimctions are analyzed in terms of density of states as sum of 
partial densities, strength functions, NPC and occupancies of single particle orbits and correlations between them. 
These results form Section IV of the paper. Section V discusses the structure of the Hamiltonian matrices and 
the modifications to be incorporated in EG0E(l+2) for a more quantitative description of chaos and localization in 
lanthanide atoms. For completeness. Section II gives the method of atomic structure calculations and Section III gives 
the definition and a brief discussion of some of the basic results for EG0E(l-|-2). Finally Section V gives conclusions 
and future outlook. 



II. METHOD OF ATOMIC CALCULATIONS 

A. Multi-Configuration Dirac-Fock method 

The earlier works on lanthanide atoms [20] have shown that the multiconfiguration Dirac-Fock ( MCDF) method 
is suitable for studying the structure and calculating properties of these atoms. The method is the relativistic 
equivalent of multiconfiguration Hartree-Fock (MCHF) [21] used extensively in non-relativistic atomic calculations. 
For a detailed description of MCDF see ref [22] and several groups have implemented it to study atom and molecules 
[23]. For completeness and continuity important features of the method are described in this section. The essence of 
MCDF is the calculation of orbitals in a mean field arising from a linear combination of configurations. Then, the 
single particle states are the eigenstates of variationally calculated single-electron Schrodinger equations and these are 
represented as \nK), where n is the principal quantum number. The quantum number n = ^{l ± s) + 1/2 is similar 
to the total angular momentum quantum number j = I ± s, however unlike j it is unique for each orbital symmetry. 
The single electron wavefunction is the two-component Dirac spinor 

..MH.K)^i(j:S::^i«,) (d 

where PnKir) and QnK,{T) are large and small component radial functions respectively, and XKm{(), 4>) ^^nd x-nmiO, 4>) 
are the spinor spherical harmonics in the Isj coupling scheme. The many electron basis functions of the method 
are the configuration state functions (CSF's) [21]. Each CSF is represented as \jPJM), where P, J and M are the 
parity, total angular momentum and magnetic quantum numbers, and 7 is a quantum number to identify each CSF 
uniquely. The CSFs are constructed in two steps: first couple the identical electrons in each sub-shell riiKi to give Xi, 
and second couple to obtain J. The construction of the CSFs are better described in second quantized notations, 
let al^j^ and anjm be creation and annihilation operators of the single electron states. Then, the al^j^ operators form 
a complete set of spherical tensor of rank j. However, for the annihilation operators 

O'jm — ( 1)'' 0'j,—m (2) 

forms a complete set and not anjm- The operators satisfy the modified anticommutation relation 

\ano X alj'V^ + i-^y^''~'[alj' x anj^^ = ^/2jTl6jj>5joSMoSnn', (3) 

where [a,y x a^ijiy'^ represents tensor coupling of the operators to rank J and component M. In absence of external 
magnetic field, like in the present calculations, the CSFs are degenerate in M and each CSF can be identified without 
M as \^PJ). A CSF I7PJ) having p sub-shells is created from the vacuum state 

I7PJ) = [(at)«^^.,X,[(at_l)«-^^.,_lX,_l[...[(a^)«^^/l,Xl]■^^]...]■^''-]■^lO) (4) 

where qi is the number of electrons in the i^^ sub-shell and Ui is the seniority quantum number, which identify identical 
sub-shell total angular momentum Xi uniquely. It is also possible to use an ordering index instead of u for details; 
see ref [24]. 
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An appropriate rclativistic Hamiltonian to describe high Z atoms hke lanthanides is Dirac-Coulomb Hamiltonian 
H^'^, which includes only the electrostatic interactions. For N electron atom 

i=l i=lj=i+ll'^» 

where oti and [3i arc Dirac matrices, pi is the electron momentum, Z{ri) is the nuclear-charge at r and N is the number 
of electrons. The first two terms are single electron Dirac Hamiltonian, the third and last terms are electron- nucleus 
and electron-electron Coulomb interactions respectively. The H^'~' is diagonal in the total angular momentum J states 
\TPJ), the atomic state functions (ASF). Here, F is a quantum number to identify each ASF uniquely and parity P 
is a good quantum number as is invariant under parity transformation. The ASFs are linear combination of 

CSFs \^rPJ) 

\TPJ)=Y,Crr\lrPJ), (6) 
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The ASFs satisfy the Schodinger equation 

H^^\TPJ) =Er\TPJ) (7) 

where is the eigenvalue. 

The starting point of the MCDF method is the variational optimization of an energy functional defined in terms of 

H^'^ and Lagrange multipliers with respect to one or more ASFs. The parameters of variational optimization are the 
coefficients c-y^ and orbitals ipnK{f)- In general, the energy functional is extremized with respect to the ground state. 
Another class of MCDF calculations is the extended optimal level (EOL), where the energy functional is extremized 
using a set of ASFs. We use the later method in our calculations and the advantages over the general method is 
explained in the next subsection. The energy functional of MCDF-EOL calculation is 

ric n^^, ric n,,; — 1 

W^DC ^ g ^^^^DC + ^ g drrq{a)ea + ^ E 5abeabN{ah), (8) 

r,s a=l r=l a=l 6=a+l 

where drs and drr are the weight factors, H^p are the matrix element of H^'^^ between the CSFs \^rPJ) and \^sPJ) 
and riw is the number of the orbitals. The quantity N{ab) is the overlap integral between the o**^ and 6*^ orbitals 
and Ca and eab are the Lagrange multipliers to enforce orthonormality between orbitals of the same symmetry but 
different principal quantum numbers. The weight factors can be chosen in several ways, in EOL calculation 

% 

where nj, is number of ASFs. Thus, atomic structure and property calculations are carried out using a hierarchy of 
eigenstates: orbitals, configuration state functions and atomic state functions. 



B. Orbital calculation and configuration interaction 

The calculation of an appropriate set of orbitals {V'(^)} is crucial in atomic structure and properties calculations. 
In the present work we calculate the orbitals using the MCDF method described in the previous subsection. This 
method can include strong configuration mixing, however the orbitals generated are state specific, this is undesirable 
for statistical studies involving several excited states. To strike a balance between these two features of MCDF method, 
orbitals are generated within selected CSF space and the energy functional is optimized with respect to several states 
using MCDF-EOL method. The orbitals generated are less state specific as these are calculated in the potential of 
several states self-consistently. The orbital set of each atom are generated in a series of calculations described in this 
section. 

The calculations are relativistic, however for compact notation we define configurations in non-relativistic notations. 
The Xe like core is considered as reference state |$o) and in non-relativistic notations 



l*o) = [(4/[(4J'[---[(aL)']]---]]|o). 



(10) 
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TABLE I: The number of relativistic CSFs iVcsF arising from each of the non-relativistic configurations. The range of the 
diagonal elements or €™'° and e™^ are also given, the ek values are defined relative to the lowest ek in the CSF space. In the 
table, the configurations are defined relative to the ground configuration 4/'"6s^, where m is 4, 5 and 6 for Nd, Pm and Sm 
respectively. 



Nd Pm Sm 



Config 


-^CSF 


mm max 


-^VcSF 


mm max 


A^CSF 


mm max 


4/ 6s 


19 


0.0 -0.25 


29 


0.0 -0.28 


46 


0.0 -0.38 


4/'"6s5d 


288 


0.31-0.64 


521 


0.09-0.53 


748 


0.03-0.58 




40 


0.39-0.69 


92 


0.33-0.65 


172 


0.39-0.85 


4/™-5d2 


608 


0.63-1.08 


1090 


0.24-0.69 


1579 


0.18-0.86 


4/'"-i6s5d6p 


684 


0.67-1.04 


1619 


0.43-0.87 


3014 


0.44-1.00 


4/'"6p2 


236 


0.97-1.30 


409 


0.83-1.25 


602 


0.87-1.38 


4/"'+i6p 


172 


1.11-1.50 


241 


1.02-1.40 


276 


1.09-1.61 




31 


1.14-1.30 


91 


1.09-1.42 


234 


1.19-1.58 


4jm-2g^25^2 


78 


1.38-1.64 


240 


1.01-1.34 


608 


1.03-1.56 


4ym+2 


46 


1.56-1.90 


46 


1.61-1.97 


46 


1.74-2.10 



where |0) is the vacuum state. The intermediate angular momenta Ji are zero as all the sub-shells are completely 
filled. The orbitals of the ground state configuration [Xe]6s^4/'"5d" are then generated treating 6s, 4/ and 5d orbitals 
as active valence shells. The calculation is in the CSF manifold of all single and double replacements among the active 
valence sub-shells of the ground configuration. The CSFs considered have the form 

\liPiJi) = Kayo's Z/5d^5d[(4j'%'^6«X6«[(4^)'-%Z/4/X4/]-^"]-^-]-^^|$0), (H) 

where "Pi , qi and Vi are the sub-shell occupancies and satisfy the condition -I- -I- = 2 -\- m -\- n. Further, the 
conditions 'm>ri>'m — 2,2>qi>0 and n > pi > n — 2 arc imposed to select single and double replacement CSFs 
within the ground configuration. Then, a MCDF-EOL calculation of the the lowest multiplet levels is carried out. 
Valence orbital 6p is generated in another calculation with frozen core approximation, where the orbitals of the ground 
configuration generated in the previous calculation arc held fixed. The CSF space of the calculation consists of single 
excitation from the ground configuration to 6p. The orbital 5d is generated in a similar calculation if it is unoccupied ( 
n = 0) in the ground configuration. Thus, in relativistic notations, the orbital set consists of (1 — 6)si/2, (2 — 6)pi/2.3/2) 
(3 — 5)^3/2,5/2 and 4/5/2^7/2- Out of these 6si/2, 4/5/2^7/2, 6^1/2,3/2 and 5^3/2,5/2 are valence shells, then the number 
of valence shells Ny and core shells Nc are 7 and 17 respectively. All the orbitals are made spectroscopic, that is, each 
orbital has n — I — 1 ( where n and I are the principal and orbital angular momentum quantum numbers) number 
of nodes. The other possibility is calculating the orbitals without imposing the constraint on the number of nodes. 
Orbitals of this type are referred to as correlation orbitals. These are usually contracted, state specific and represent 
correlation effects of specific states very well. For statistical properties calculations, the spectroscopic orbitals are 
appropriate. The MCDF calculations to generate the orbitals include important intra valence correlation effects but it 
is incomplete within the orbital basis considered. A configuration interaction (CI) calculation, which is diagonalization 
of H^'-^ matrix, of all the possible CSFs can capture the remaining correlation effects. However, such a calculation is 
impossible as the size of the CSF space increases exponentially with the number of active valence electrons and Ny. 
In our calculations the CSF space consists of single and double excitations from the ground configuration 

\jM) = [{alr%iyepXep[{4,r%i^,dX54^^^^^^^^^ (12) 

where Pi + Qi + ri + Si = 2 + m + n, and impose the conditions m>Si>m — 2, 2>ri>0, n>qi>n — 2 and 
2 > Pi. Thus, all the CSFs considered are connected to the ground configuration [Xe]6s^4/'"5d". The CSF space has 
ten non-relativistic configurations, among these 4f'^5(P and 4/™~^6s5(i6p together contribute more than half of the 
CSFs. The maximum number of CSFs arises from 4/™^^6s5(i6p, which has the largest number of open shells. The 
number of CSFs arising from each non-relativistic configuration and range of {jPJ\H^'^\'yPJ) are given in Table I. 

Consider A^^ and as the number of active valence electrons and single particle states respectively. Then, Na is 
6, 7 and 8 for Nd, Pm and Sm respectively and A^^ = 32 since the valence space consists of 6si/2, 6^1/2^3/2, 5^3/2,5/2 
and 4/5/2,7/2 ( each orbital contribute 2j + 1 number of single particle states ). In addition, define the number of 

unfilled single particle states as Ng = Na — Na. Then, the number of possible determinantal states is ( "° ) and the 
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two-electron Coulomb interaction couples each Slater determinant to 

K = l + NaNo + Na{Na - l){No){No - l)/4 (13) 

determinants which include itself and others with shell occupancies different by one and two. In a single determinant 

approximation of the ground state, K is also the number of determinants in the manifold which inchidcs the single 
and double excitations. This is the determinantal equivalent of the CSF manifold chosen in our calculations. For 
further analysis, consider these K determinants as the manifold chosen for calculations. As mentioned in the context 
of CSFs, the representation of two-electron Coulomb interaction is incomplete in such a determinantal manifold. This 
is because, the Coulomb interaction couples pairs of determinants which differ by one or two occupancies. Hence, 
for the singly excited determinants, the zeroth ordc;r dc;scription of Coulomb interaction is incomplete without the 
triply excited ones. The quadruply excited are similarly needed for the doubly excited determinants. Among the K 
determinants, the fraction coupled to a singly excited one is 

^ K-{N,- l)(iVa - 2)(iV„ - 1)(A^, - 2)/4 ^^^^ 
~ when N„ » N^. 

The last term in the numerator is the number of triply excited states coupled to a singly excited determinant but 
absent in the manifold. Similarly, the fraction of the determinants coupled to double excited determinant is 

^ K - [{Ng - 2){No - 2) + {Ng - 2){N, - 3)(iVo - 2){No - 3)/4] 

= (15) 

~ when N^:^Na. 

That is, each of the states considered are directly coupled to ~ l/iV^ of the total number. However, higher order two- 
electron Coulomb interaction connects all the determinants within the same J manifold. For Ng ^ 3, the Kg and Kj^ 
scaling show the number of determinants mixed at first order is close to the GOE predicted value of number of principal 
components (NPC). Similar scalings apply to the CSFs as these are linear combinations of Slater determinants. This 
implies that a CSF space consisting of single and double excitations alone can provide a good representation of the 
eigenvalues and eigenfunction properties. 



III. EIGENFUNCTION PROPERTIES IN EGOE(l+2) 



This section gives a brief introduction to EG0E(l-|-2) and a summary, for later use, of the results known for this 
random matrix ensemble; see [11] for a review on this subject and also [18]. 



A. Definition of EGOE(l+2) 



Let us start with two-body embedded GOE, i.e. EG0E(2) or TBRE defined for spinless fermion systems. The 
EG0E(2) for m {m > 2) fermion with the particles distributed say in N single particle states | i^i) , i = 1, 2, . . . , iV) is 
generated by defining the Hamiltonian H, which is 2-body, to be GOE in the 2-particle space and then propagating 
it to the TO-particle spaces by using the geometry (direct product structure) of the m-particle spaces. Operator form 
for a 2-body Hamiltonian H=V(2) is defined by, 

^(2) = J2 ^""^ I ^(2) I yj)g al^ a,, a,, ; (16) 

{vk VI I V{2) I Vj Vi)^ = - {vk n I V{2) I Vi Vj)^ , 
{uk vx I V{2) I Vi Vj)^ = {ui uj I Vi2) I n)g 

Now, the Hamiltonian H matrix in m-particle spaces, in the occupation number basis (occupation numbers will be 

or 1 and the basis is generated by distributing the m particles in all possible ways in the single particle states 
fr), is defined in terms of the two-body matrix elements {v^ vi \ V{2.) \ Vi Vj) g (note that the subscript 'a' stands for 
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antisymmetrized two-particle states) and the non-zero matrix elements are of three types, 

(l/l!/2 • • • I ^(2) I ■ ■ ■ Vm)a = ^^^^3 ' I ^i^i") a 

Vi < Vj < I'm 

■■■Vm I ^(2) I l'lV2 ■ ■ ■ l^m)a = ^ i'^p'^i I ^(2) I t^lt^i) a 
{VpVqV^ ■■■Vjn \ V {2) \ ViV2Vz ■ ■ ■ Vm) ^ = {I'p'^g \ V{2) \ l'll'2) ^ I (17) 

all other (•••|_ff|---)^=0 due to the two-body selection rules. 

The EG0E(2) is defined by the above Eqs. and a GOE representation for V(2), 

I ^(2) I Vi Vj)^ are independent Gaussian random variables 

{vk VI I V{2) I V, v^)^ - , 

\{v^ VI I V{2) I I.,.) J' = v\l+5^i^)^(kl)) (18) 

where the bar denotes ensemble average and v"^ is a constant. Note that d{m) = ( " ) is the Hamiltonian matrix 

dimension and the number of independent two-body matrix elements is [d{2){d{2) + l)]/2. 

Hamiltonian for realistic systems such as atoms consists of a mean-field one-body (defined by a finite set of single 
particle states) plus a complexity generating two-body interaction. Then the appropriate random matrix ensemble, 
for many purposes, is EG0E(l-|-2) defined by 

{H} = h{l) + X{V{2)} (19) 

where { } denotes an ensemble. In (19), the mean-field one-body Hamiltonian h{\) = J2i^i''^i ^ fixed one-body 
operator defined by the single particle energies with average spacing A (note that is the number operator for 
the single particle state In general one can choose Ci to form an ensemble; see [14] for examples. The {y(2)} 

in (19) is EG0E(2) with = 1 in Eq. (18) and A is the strength of the two-body interaction (in units of A). 
Thus, EG0E(l+2) is defined by the four parameters (m, TV, A, A) and without loss of generality we choose A = 1. 
Construction of EG0E(l-|-2) follows from Eqs. (16), (17) and (18) by just adding to the diagonal matrix elements 
(see the first equality in (17)) the term e,^. where Vi are the occupied single particle states for the given basis state. 



B. Basic properties of EGOE(l-|-2) 

Most significant aspect of EG0E(l-|-2) is that as A changes, in terms of state density, level fluctuations, strength 
functions and entropy, the ensemble admits three chaos markers as described in Fig. 1 and reviewed in [18]. Firstly, 
it is well known, via the EG0E(2) results in [10, 25] and the fact that general /i(l)'s produce Gaussian densities, 
that the state densities p^'^{E) = {d{H — E))"^ take Gaussian form, for large enough m, for all A values (often the 
superscripts {H,m) are dropped), 

p"'"^{E) = \ exp (20a) 



E={E-eH{m))/aH{m) (20b) 

In (20), eH{rn) = (H)™ is the spectrum centroid and similarly aij(m) is the spectrum width. In practice there will 
be deviations from the Gaussian form and they are taken into account by using Edgworth expansion [26] in terms 
of the skewness (71) and excess (72) parameters. With {E)dE = r]{E)dE, the r]g{E) with lower order Edgeworth 
corrections is, 

Ved{E) = 

Vq{E) |i + ^-^He^{E) + f^He^{E) + ^^He,{E)^ ^^^^ 
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where Her{E) are Hermite polynomials. With A increasing, there is a chaos marker Ac such that for A > Ac the level 
fluctuations follow GOE, i.e. Ac marks the transition in the nearest neighbor spacing distribution from Poisson to 
Wigner form. Parametric dependence of Ac is Ac oc l/m?N [27]. Now let us define strength functions. 

To maintain consistency with the notations used in EG0E(l+2) literature, from hereafter, we represent CSFs 

Given the mean-field h{l) basis states 
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as |A;) and ASFs as \E). Then, Eq. 6 assumes the form \E) — '^^cf\k). 
1^) = YliE 1-^)' strength functions (one for each k) Fk{E) = 



E,I3 



{dp"{E)). As A 

increases further from Ac, the strength functions change from Breit-Wigner (BW) to Gaussian form and the transition 
point is dented by Xp; see Fig. 1. The BW and Gaussian (denoted by Q) forms of Fk{E) are, 
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Fk-.BwiE) = — 



2n {E 



Ekr + rl/4' 



(22a) 



{E - Ek f 



F^..g(E) = -^=^ exp , (22b) 



where E^ = {k\H\k). With p = J^"^ Fk{E)dE, the spreading width = £3/4 — S1/4. Similarly the variance of Fk 

is al = (k \ H'^ \ k) — {{k \ H \ fc))^. For Ac < A < Af (this is called BW domain) they are close to BW form and for 

X > Xp (this is called Gaussian domain) they approach Gaussian form. In fact the BW form starts in a region below 
Ac (as shown in Fig. 1, there is a Aq such that below Aq, the strength functions are close to 5- function form and for 
A > Ao there is onset of BW form) but the fluctuations are close to Poisson for A < Ac. Arguments based on BW 
spreading widths give Af oc 1/y/m [5, 28] and this result will be used later. Unlike the Poisson to GOE transition 
in level fluctuations, which is common for Hamiltonians with regular and irregular (chaos generating) parts, the BW 
to Gaussian transition in strength functions is a unique signature for the operation of EG0E(l-|-2) in finite quantum 
systems. For the BW to Gaussian transition, a good interpolating function has been constructed recently [18], 

Fk..BW-g{E : a, f3) dE = 

r(a) dE (23) 

V^r(a-i) HE-Ekr+aP)- 

In Fk.Bw~g{E : a,/?), /? supplies the scale while a defines the shape (hence a is the significant parameter) as it gives 
BW for a = 1 and Gaussian for a ^ 00. The variance of Fi^-Bw-g is <y'^ {Fk-BW -g) = '^1 — ^(^/ (^a — 3), a > 3/2. It is 
important to stress that in general (a, /?) change with k although for EG0E(l-|-2) they are nearly constant except for 

with Ek very far from their centroid. Fig. 2 shows that the interpolating function (23) fits very well the numerical 
EG0E(l-|-2) results for strength functions; for more detailed discussion see [18]. 

One important measure of the complexity of eigenstates of interacting systems is the participation ratio PR 
(denoted by ^2), £,2{E) = {X]fe l^'fe'l^}""'^ where the subscript '2' for ^ denotes that it is the second Rcnyi entropy [29]; 
in literature PR is also often called number of principal components (NPC) as it gives the number of basis states that 
make up the eigenstate with energy E. Following Ref. [30] one can write ^2 in terms of the strength functions, 

{6(i?)/^2^°^}" =7^7— 2 / dE,p^{E,)[F,{E)]\ ^f°^ = d/3 (24) 
[P (E)] J-00 

In writing (24), it is assumed that the level and strength fluctuations are of GOE type and hence it is valid only for 

A > Ac- Further, p^{Ek) is the density of Ek energies and it is generated by an effective one-body Hamiltonian h [30]. 
Also, in general p^{Ek) is also a Gaussian for EG0E(l-|-2). Assuming that a and /3 are fc-independent, £,2{E) in the 
BW to Gaussian transition domain can be evaluated by substituting Fk^Bw-g for Fk{E) in (24). At the spectrum 
center the integral in (24) can be evaluated and it gives [18], 



UE = 0)/^c 



GOE _ 



I Y (2a - 3) r2(a - i) VC'(1 - C^) (25) 
(1,1 -2a, i^^^-^)}"^ 



where U{ ) is hypergeometric-U function [31]. In (25), the correlation coefficient C = yl — ffe/cl^ where 

i-/'^) ^i^j ^ij '^H — (i-/"^) J2i{Ei — e)^ are the variance of the off-diagonal Hamiltonian matrix elements 



'k 
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and eigenvalues respectively. However, in the Gaussian domain, with F}. [E) being a Gaussian, gives for any E, 



um^'"'' = «^p-rTc^ (26) 

Instead of NPC, it is possible to use the closely related information entropy S'^'f^iE); S'^^^^iE) = 
\^k\ ^'^sl^fe'l ■ interesting recent observation is that [32] as we increase A much beyond \p, there is 
a chaos marker Af around which different definitions of entropy (for example 3™-^°, thermodynamic entropy defined 
via p^{E), single particle entropy defined via occupation numbers), temperature etc. will coincide and also strength 
functions in h{l) and V{2) basis will coincide. Thus A At region is called the thermodynamic region and [13] gives 
first application of this marker. 

It is important to point out that the cigcnstates of atoms carry good angular momentum (J) (in some situa- 
tions even good L, S) and therefore in principle one should consider EG0E(l+2) preserving J-symmetry (called 
EG0E(l+2)-J in [11]). Theory for EG0E(l+2)-J is not available (see however [33, 34] for some first attempts), and 
just as in the nuclear shell model studies [11, 35], it is assumed that the forms for p^{E), Fk{E), (NPC)^; etc. derived 
using EG0E(l-|-2) for spinless fermion systems, extend to J spaces. This can be considered to be an aspect of chaos 
in atoms and this is the basis for statistical atomic spectroscopy as being developed by Flambaum; see [3, 8, 36] and 
references therein. 



IV. RESULTS FOR p"iE), Fk{E), ^2(£) AND (n„>^ 

In the present study, we select the J = 4+ CSF manifold for Nd and Sm atoms. The parity is chosen equal to the 
parities of the ground configurations 4/^6s^ and 4/''6s^ respectively and the total angular momentum J is the same 
as in the Ce I work [2]. For Pm, which has an odd number (i.e. 7) of active valence electrons we choose J — 9/2~ 
CSFs. The parity is again equal to the ground configuration parity 4/^6s^. It is to be reiterated that, the orbitals 
of the present study are different from our earlier work [17]. As mentioned in Section II B, the orbitals in this work 
are calculated in a multiconfiguration potential, whereas in [17], the orbitals are calculated in a single configuration 
potential. The manifolds of the Nd, Pm and Sm have 2200, 4375 and 7325 CSFs respectively. The CSFs are generated 
in the sequence of sub-shell occupations but later, these are energy ordered. The CI calculation is carried out within 
the whole space for Nd and Pm. For Sm, as mentioned in [17], the density of states p{E) is bimodal. We select the 
first 6300 CSFs which contribute to the first peak for the CI calculation. The atomic state \E) of the Nd, Pm, and 
Sm obtained from the CI calculations have energy range of 2.20, 2.27 and 1.60 hartree respectively. 



A. Density of states p^ (E) 

By definition, density of state p{E) = {5{H — E)) is a sequence of 5-functions. For obtaining the smoothed (with 
respect to energy E) density of states and for comparison with the EG0E(l-|-2) predictions of p^-"^{E), we construct 
the binned density of states. In the our calculations, the energy range of \E) of each atom is divided into sixty bins 
of equal sizes. The density of states is then binned and normalized, these are shown in Fig. 3. The Gaussian and 
Edgeworth corrected Gaussian calculated from the centroid e, variance cr^, skewness 71 and excess 72 of p{E) are also 
shown in the figure ( e, cr^, 71 and 72 values are given in the figure). For SmI Gaussian is in good agreement with 
the p{E) except at and around the centroid. At the centroid the Gaussian is ~ 13% less than the actual p{E). In 
general all p{E) are multimodal in structure and it is most prominent in Nd. The p{E) for Nd has a pronounced peak 
at E ^ —1.5 and fluctuations above the mean for E between 1.0 and 2.5. Similarly for Pm, p{E) has a pronounced 
peak at ii^ ~ 1.4. As discussed ahead, the multimodal structures correlate with similar structures of the number of 
principal components and observed deviations of the strength function from the dependent variation. Finally for 
SmI, there is a peak a.t E ^ 2. 

The multimodal form of p^'"^{E) can be understood by decomposing it into a sum of partial densities. For each 
sub-shell occupancy m one can define a partial density p^'"^{E) with centroids Ec{fh) and a'^{m) defined by, 

E,{fh) = {Hf , a^{m) = {{H - E,{m)ff . (27) 
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Note that cr^(m) is sum of partial variances a^{m — > m') where 

CT2(m) = ^ (T^(m m!) , 

rh' 

a'^im m') = {d{m)y^ ^ \{ma\H\m' p)f for mj^m' 



a, 13 



(28) 



a2(m ^ m) = |rf(m)-i ^ |(m a | | m - {Ec{m)f 

Usually the (to fh) is called internal variance and others, partial variances [37] . In terms of the partial densities 
the total density of states is, 

p"'"'{E) = [d{m)]-^ ^p"'^{E) d{m) (29) 

where to is the number of active valence electrons, to represent the possible sub-shell occupancies, ^(to) and ^(to) are 
the number of CSFs in the CI calculation and number of CSFs in the to sub-shell occupancy manifold respectively. 
That is, p^'"^(E) is the sum of partial densities p^'™-{E) from the sub-shell occupations defined in Eq. (12). It is 
important to stress that Eq. (29) is exact. As evident from Eq. (12), each of the possible sub-shell occupations 
have different angular momentum coupling sequence. It is easily seen that p^'"^{E) takes multimodal shape when 
the partial densities are well separated with weak mixing between them (i.e. partial variance is much smaller than 
the internal variance). The range of ct (i.e. diagonal elements of the Hamiltonian matrix ) of each non- 

relativistic configuration given in Table I, which is a representation of the partial density spread, is plotted for Nd 
in Fig. 4a. At ~ 0.62 hartree, it shows small overlap between the Af^Qsbd and next in energy 4/'*5(i^ and 4/^6s5d6p 
configurations. The separation is also visible in the Hamiltonian structure discussed ahead [Fig. 9 ahead shows two 
distinct blocks centered around ~ 0.5 hartree and ^ 0.8 hartree respectively]. In addition, the Coulomb mixing 
(4/''6s5d|l/ri2|4/^5d^) and (4/^6s5d]l/ri2|4/^6s5(i6j?) are weak as the configurations differ by a single occupancy. 
The separation between the configurations persists in \E) and the contributions from 4/^6s5rf configuration manifests 
as the smaller peak of the bimodal p^'"^{E). The energy range of the remaining configurations show large overlaps 
and contribute to the main peak of the Nd p^'"^{E). The 4/^6s^6y» configuration which lie between 4/^6s5d and 
4/^5d^ energetically, is single replacement from 4f'^5d'^ and doesn't contribute to smoothing the p^'"^{E). Fig. 4b 
shows the result of Eq. (29) with partial densities represented by Gaussians for Nd. The agreement between theory 
and calculations is excellent. This demonstrates that it is possible to construct the total density of states as sums 
of Gaussian partial densities (even though a single Gaussian representation as shown in Fig. 3 deviates strongly 
from the calculated results). This result is significant as it is possible to calculate the centroids and variances (also 
partial variances) without recourse to H matrix construction [37]. Similar to the Nd, the bimodal structure of the 
Pm's p^'™'{E) is a consequence of the small overlap between the partial densities of the 4/^6.s5d6p aiid Af^Qp^ 
configurations. However, unlike Nd, these configurations are in the higher end of the energy range. In contrast to Nd 
and Pm, the Sm atom has partial densities with large overlap and the p^''^{E) is without prominent secondary peaks. 
Now we will turn to strength functions which are nothing but the partial densities with to's divided into individual 
basis states. 



B. Strength functions Fk{E) 

The Fi.{E) of the individual CSFs \k) exhibit large fluctuations and significant variation is observed between Fk{E) 
of neighboring \k). This is evident from the selected \k) presented in our Sm I work [17]. For statistical description 
we consider the representative strength function Fk{E) which is the average of individual Fk{E) within the range of 
energy efe ± Ae. In the present work, we have chosen efe ± .25(Tfe as the range of averaging to calculate Fk{E) and the 
centroids of the individual Fi.{E) are aligned while averaging. Fig. 5 shows F^^E) of Nd I, Pm I and SmI calculated 
at Cfe = —1.5, —1.0, —0.5, and 0. For all the three atoms, the interpolating function of Eq. (23) fits the Fk{E) very 
well; the f3 parameter in Eq. (23) has been eliminated using the calculated cTfe. The best fit value of a for Nd I, Pm 
I and SmI for ?fe = are 5, 7.5 and 11 respectively, thus these are in the Breit-Wigner to Gaussian transiti on regio n 
with SmI closest to Gaussian. This result is already reported in [18] (it should be noted that in this paper Fk{E) is 
constructed by averaging Fk{E) over 3% of the basis states around e'k = and this is different from the procedure 
used in Fig. 5 and therefore there is slight difference in the a values extracted). Moving away from 'ck = 0, firstly it 
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is seen that the parameter a in Eq. (23) is k dependent and it changes from a value close to that of BW to a value 
giving close to Gaussian form as is increasing from -1.5. For Nd, Pm and Sm a changes from 2 to 5.2, 2.5 to 7.5 
and 2.5 to 10.6 respectively. The variation in a seen in Fig. 5 is understood from the fact that near the Fermi surface 
(i.e. for Cfc ~ —1.5 to —1) the levels are well separated and hence the mixing is weak giving BW form and as we go 
towards the center the mixing is strong giving close to Gaussian (or between BW and Gaussian) form. It should be 
pointed out that EG0E(l-|-2) ensemble do not produce large changes in a as /c changes in Fk{E). However, in all 
the calculations it is seen that ak is essentially constant as expected from EG0E(l+2). Clearly, the Fk{E) analysis 
shows that modifications of EG0E(l-|-2) are needed. 



The inverse participation ratio ^2 {E) of an eigenstate is the effective number of basis functions contributing to it. It 
provides a measure for the presence of chaos in the system. For a GOE with dimension d, the NPC is d/3 independent 
of energy while for EG0E(l+2) there is strong energy dependence (it is Gaussian in the Gaussian domain; see Eq. 
(26)). The ^2iE) of Nd, Pm and Sm for all \E) are calculated and shown in Fig. 6a-c. The ^2iE) of Nd and Pm 
are multimodal and localized states (i.e. states with NPC much smaller than GOE or EG0E(l-|-2) predicted values) 
are prominent in all the atoms. We should add that localized states are also observed in the previous studies on 
Ce I [2] and Pr I [4]. In addition, more recently, localization quite similar to those shown Fig. 6 are also seen in 
J7^^+ calculations [38]. Despite localization, Nd and Pm have Breit-Wigner like envelopes, whereas for Sm it is like 
Gaussian. Nd and Pm have Gaussian like secondary peaks around ~ — 1.5 and ~ 1.4 respectively and it is to be noted 
that, these are the locations where p^'"^{E) has local peaks. For Sm, like p^'"^{E), it is unimodal. At the centroid 
of the envelopes, the value of (,2{E) are ^ 0.39, ^ 0.42 and ^ 0.5 for Nd, Pm and Sm respectively. These should be 
compared with the EG0E(l+2) values given by Eq. (25) (for the Nd, Pm and Sm atoms, is 0.88, 0.86 ans 0.83 
respectively) and they are 0.44, 0.48 and 0.55 respectively. Here it is assumed that a value at the = can be 
used for all e^'s. However this is not a good approximation as seen from Fig. 6. In a better calculation we use Eq. 
(24), where the integral is divided into seven segments each with spread .5 in the range -1.75 to 1.75. Within each 
segment Fk{E) is constructed with a taken for the mid e^. but incorporating exact Fk{E) centroids and widths. This 
calculation gives overall good description for Sm while for Pm there are deviations for > 0. This is not surprising 
because there is a secondary peak for E > 0, which cannot be accounted in the EG0E(l-|-2) model adopted. For Nd 
we have not shown the result as the a for ik = —1.0 is undetermined. Assuming an interpolating value of a, it is seen 
that the calculations describe reasonably well ^2(-E') for < 0. 

In Fig. 6d-f ^2 {E) is binned to remove local fluctuations and the resulting histograms are shown. They give at the 
center the (,2{E) to be 0.27, 0.33 and 0.48 for Nd, Pm and Sm respectively. Finally, though the number of localized 
states appear numerous in the plots, they are in fact not that many. This is seen clearly from Fig. 6g-i where the 
relative density of localized states is plotted as a histogram. We classify a state as localized if it has ^2 (E) less than 
or equal to 20% of the binned value in corresponding energy range. As shown in Fig. 6g-i the number of localized 
states in each energy bin pioc{E) is small for all the atoms studied. This indicates the difference between the envelope 
and binned ^2iE) is due to high frequency fluctuations and not due to purely localized states. In order to understand 
more about localized states an analysis is carried out in terms of occupancies of single electron orbits and we will turn 
to this now. 



In the statistical description of a many-body quantum system, properties of the system are often expressed as 
functions of the sub-shell occupation numbers [8] . In this section, we use sub-shell occupation numbers to study the 
nature of the localized states. By definition the occupation number of the sub-shell a in |A;) is 



C. Number of principal components ^2iE) and localized states 



D. Occupation number {ua}^ and the localized states 





(31) 



k 



It is evident from (,2{E) calculations that, there are several localized states as shown in Fig. 6g-i. Each of these have 
^ioc(-E) <C ^2{E) and few can represent these states. A possible origin of localized states is strong Coulomb mixing 
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between few |fc) of similar sub-shell occupancies, then the two-electron Coulomb integral is between same orbitals 
and is large { K — multipole is allowed). This is indeed the case and the signature is the anticorrelation in the 
trend of £,2{E) and occupation number of selected sub-shells shown in Fig. 7. The figures show £,2{E) and (ua)^ of 
few \E) around centroid energy. In Fig. 7a, there is anticorrelation between the ^2{E) and {n^d^fz)^ + (^5^3/2)^ of 
Nd, similar trend is observed for Pm in Fig. 7b i.e. the occupancies are large when NPC is small. This shows the 
localized states of Nd and Pm within the spectral range considered arise from strong mixing between of high hd 
occupancy. Whereas, the localized states of Sm as shown in Fig. 7c arise from the strong mixing between \k) of high 
6s or 6p occupancies. 

To compare and contrast the chaotic and localized states, |C^p of pair of neighboring chaotic and localized states 
of Nd, Pm and Sm are shown in Fig. 8. The figures clearly shows that the localized states ( energy states, T = 800, 
r = 2148 and F = 2885 of Nd, Pm and Sm respectively) are represented by few \k), whereas the neighboring chaotic 
states ( F = 801, F = 2149 and F = 2886 of Nd, Pm and Sm respectively) has contributions from several |fc). The 
|C^P of the chaotic states exhibit a systematic growth and decay, in contrast localized states show no systematic 
trends. To show the anticorrelation between {ua) of selected sub-shells and S,2{E), the contributions from each basis 
|C|'p(na)fc are plotted. For the localized states the plots are for a that has the highest value. For the chaotic states, 
all the sub-shells show similar trend and approximately an order of magnitude less ( see the scales in Fig. 8) than 
that of the localized states. 



V. HAMILTONIAN MATRIX STRUCTURE 

Understanding the structure of the Hamiltonian matrix is crucial for developing an appropriate random matrix 
model for the Lanthanide atoms. The diagonal Hamiltonian matrix elements H^j^ will have contributions from both 
the one and two particle Hamiltonian terms. However, the off diagonal terms have non-zero contributions 

from the Hamiltonian terms depending on the relative excitation between \k) and \k'). If the two CSFs are singly 
excited with respect to each other, then one and two particle terms contributes, whereas only the two particle term 
contributes when the two CSFs are doubly excited with respect to each other. The off diagonal matrix elements are 
zero for triple or higher excitations. As a consequence, the Hamiltonian matrix is sparse. In our present calculations, 
several of the CSFs which are doubly excited with respect to ground configuration are quadruply excited with respect 
to each other. For example the configurations 4/™~^5rf"+^6s^ and 4/'"5d"'6p^ are doubly excited with respect to 
the ground configuration 4/™5(i"6s^, however they are quadruply excited with respect to each other. Hence, all the 
matrix elements between CSFs arising from 4/'"~^5(i"+^6s^ and 4/'"5d"6p^ are zero. There are other configurations 
similar to these in the CI space considered. Due to these configurations, the Hamiltonian matrix in our calculations is 
sparse, though we consider only the single and double excitations from the ground state configuration. For bringing 
out a coarse grained structure of the Hamiltonian matrix, the square of matrix elements {H^^Y (to remove the phase 
factor) is calculated and binned in terms of the diagonal elements H^k- '^^^ element of the binned Hamiltonian 
for bin size A is then 

{i'-Hkk)<^, U'-Hu)<A 

where i' = i x A and j' — j x A. The plot of binned Hamiltonian for Nd. Pm and Sm arc shown in Fig. 9. To improve 
the contrast the plots show ln{Wij) instead of the Wij. In all the three atoms: (i) band like structure is evident, 
however, there are also block structures within the band; (ii) another important feature is the prominent streaks of 
large matrix elements parallel to the diagonal, these mix two diagonal blocks and, as mentioned in [39], perturbative 
calculation can account for these distant elements; (iii) the Hamiltonian matrices of Nd and Pm have small range of 
energy with no states and these are the empty stripes in the figures. 

Presence of prominent diagonal blocks in the figures (point (i) above) is a consequence of the spread of the various 
configurations shown in Fig. 4 and Table 1. Though it is tempting to conclude that there is a BRM like structure, 
it should be noted that the various blocks retain their identity each with their own internal variance (j^(m m). 
Fig. 9 clearly shows that the partial variances cr^(m — » m') (except in some special cases) are in general much 
smaller compared to the internal variances. It should be pointed out that block matrix structure as in Fig. 9 is also 
observed very recently in nuclear shell model calculations [40, 41] and it is plausible that such structure is a feature 
of interacting particle systems. In [40, 41] it is argued that one of sources for this structure is J symmetry. The 
various blocks are much more well separated in Nd and Pm as compared to Sm and as a result Sm is much closer to 
EG0E(l-|-2) with Nd and Pm exhibiting larger localization with departures from EG0E(l+2) predictions. 

To gain better understanding of the off-diagonal matrix (i.e. (i) and (ii) above), which will mix the CSF's, let us 
consider their structure in some detail. In the present calculations the off diagonal Hamiltonian matrix elements 
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are largely generated by the inter-electron Coulomb interaction l/r\2 than the one electron terms ti where, 

U = cdi ■ Pi + c^pi - 1) - (33) 
In second quantized notation, H^'~^ can be written as 

if DC ^ V(a|i |/3)4a„ + V (75| — \aP)a\alapac, (34) 

a/3 a^"f6 

where for short we use Greek alphabets to denote orbitals and summation is without restrictions. The two-particle 
matrix element of l/ri2 can be decomposed into multipole components, 

(7(51 — |a/3) =5^G^(7,'5;a,/3)i?^(7,5;a,/3) (35) 
ri2 ^ 

where {-y,S;a, (3) and {'y,5;a, (3) are angular and radial integrals respectively, and K is the multipole. The 
angular momentum selection rules are consequence of G^{'y,5;a,(i) and the radial integral also known as Slater 
integral [42]. In total there are N^{Nc + Ny — 1)^/4 matrix elements of l/ri2 which contribute to the off-diagonal 
Hamiltonian matrix elements; A^^ and are defined in Section II. The many-particle matrix elements of l/ri2 can 
be written as, 

(7^PJM| — |7,PJM) = V Va^'-^^(7,5;a,/3)i?^(7,.5;a,/3) (36) 

al3j5 K 

where C^g '^^ (j, S; a, (3) is the angular factor which transports the single particle matrix elements to the many-particle 
space. The C"s consists of two components, first the angular factor in the single electron space and second a component 
which propagate from the single electron to many electron space. These will have complex structure when electrons 
of several open-shells are coupled and the distribution of the non-zero Hamiltonian matrix exhibit intricate patterns. 
Interestingly, nuclear shell model calculations showed that [40, 41] the angular factors generate a part of the GOE 
structure of the diagonal blocks. The radial integral is large for if = in Eq. (35) and it is an allowed multipole 
when = and Kfj = Kg- This can couple energetically well separated configurations strongly and explain (ii) 
above. Finally, from the matrix structure one can understand qualitatively the departures of Fk{E) and ^2{E) from 
EG0E(l-|-2). In this random matrix ensemble the two-particle matrix elements variance (the parameter in Eq. 
(18)) is independent of its position (i.e. independent of the indices (i,j,k,l) in Eq. (18) except that the diagonal 
matrix elements have twice the variance of the off-diagonal matrix elements). 



VI. CONCLUSIONS AND FUTURE OUTLOOK 



In this paper wavefunctions of complex lanthanide atoms Nd, Pm and Sm are analyzed in terms of strength 
functions, number of principal components, occupation numbers and also the Hamiltonian matrix structure. Exam- 
ination of Fk{E) showed that the BW form is dominant in Nd, they are more towards Gaussian in Pm and quite 
close to Gaussian in Sm for the basis states not very far from the e/c centroids. All three atoms exhibit BW form for 
Efe < —1. Thus the statistical spectroscopy developed by Flambaum et al [5-7] will be good for as long as Fk{E) with 
Cfe ^ — 1 contribute to observables of interest. In all these atoms there are localized states throughout the spectrum 
( though the density is not high) and this is clearly seen in the measure NPC. The structure of the localized states 
is closely correlated with the occupation numbers. All these and the matrix structures showed that EG0E(l-|-2) 
need to be modified. One approach is to partition the two-particle space and employ different variances for different 
block matrices and then propagate the ensemble to many particle spaces. Very few properties of these partitioned 
ensembles [p — EGOE{l -\- 2)) are known till now [11] and clearly they should be studied in more detail. The induced 
TBRE introduced recently [43] are closely related to p-EG0E(l-|-2). Another approach is consider the two-particle 
Hamiltonian as a mixture of a regular part (say determined by the lower order multipoles in Eq.(36)) and a random 
part and analyze them as a function of the strength of the random part (they are called K + aEGOE{2) in [11] 
where K is a fixed operator, and numerically some studies of such an ensemble are carried out using nuclear shell 
model recently [12]). Thus an important outcome of the present detailed analysis (presented in Sections II- V) is that 
for further progress in understanding the wavefunction structure in complex atoms it is necessary to analyse various 
modified EGOE's and this is for future. 



13 



[1] W. C. Martin, R. Zalubas, and L. Hagan, Atomic Energy Levels - The Rare-Earth Elements, Natl. Bur. Stand. Ref. Data 

Sc., Natl. Bur. Stand. (U. S.), NBS-60, (U. S., GPO, Washington, DC, 1978). 
[2] V. V. Flambaum, A. A. Gribakina, G. F. Gribakin, and M. G. Kozlov, Phys. Rev. A 50, 267 (1994). 
[3] V. V. Flambaum, A. A. Gribakina, G. F. Gribakin, and 1. V. Ponomarev, Physica D 131, 205(1999). 
[4] A. Cummings, G. O'Sullivan, and D. M. Heffernan, Journal of Phys. B 34, 3407 (2001). 
[5] V. V. Flambaum, and F. M. Izrailev, Phys. Rev. E 56, 5144 (1997). 

[6] V. V. Flambaum, A. A. Gribakina, G. F. Gribakin, and I. V. Ponomarev, Phys. Rev. E 57, 4933 (1998); A. A. Gribakina, 

G. F. Gribakin, and V. V. Flambaum, Aust. J. Phys. 52, 443 (1999). 

[7] V. V. Flambaum, G. F. Gribakin, and F. M. Izrailev, Phys. Rev. E 53, 5729 (1996); V. V. Flambaum, A. A. Gribakina, 
and G. F. Gribakin, Phys. Rev. A 58, 230 (1998). 

[8] V. V. Flambaum, A. A. Gribakina, G. F. Gribakin, and C. Harabati, Phys. Rev. A 66, 012713 (2002). 

[9] G. F. Gribakin and S. Sahoo, J. Phys. B 36, 3349 ( 2003). 
[10] T. A. Brody, ,]. Florcs, ,J. B. French, P. A. McUo, A. Pandey, and S.S.M. Wong, Rev. Mod. Phys. 53, 385 (1981). 
[11] V. K. B. Kota, Phys. Rep. 347, 223 (2001). 

[12] V. Velazquez, and A. P. Zuker, Phys. Rev. Lett. 88, 072502 (2002); V. Velazquez, J. G. Hirsch, A. Frank, and A. P. Zuker 

Phys. Rev. C 67, 034311 (2003). 
[13] V. K. B. Kota, Ann. Phys. (N.Y.) 306, 58 (2003). 

[14] Ph. Jacquod, and A. D. Stone, Phys. Rev. Lett. 84, 3938 (2000); Y. Alhassid, Ph. Jacquod, and A. Wobst, Phys. Rev. B 
61, R13357 (2000); T. Papenbrock, L. Kaplan, and G. F. Bertsch, Phys. Rev. B 65, 235120 (2002). 

[15] B. Gcorgcot, and D. L. Shepelyansky, Phys. Rev. E 62, 3504 (2000); G. P. Berman, F. Borgonovi, F. M. Izrailev, and 
V. I. Tsifrinovich, ibid. 65, 015204 (2001); V. V. Flambaum and F. M. Izrailev, ibtd. 64, 026124 (2001). 

[16] G. P. Berman, F. Borgonovi, F. M. Izrailev, and A. Smerzi Phys. Rev. Lett. 92, 030404 (2004). 

[17] D. Angom and V. K. B. Kota, Phys. Rev. A 67, 052508 (2003). 

[18] D. Angom, S. Ghosh and V. K. B. Kota, Phys. Rev. E 70, 016209 (2004). 

[19] J. Bauchc and C. Bauchc- Arnault, Comp. Phys. Rep. 12, 1 ( 1990). 

[20] D. Angom, I. Endo, A. Fukumi, M. linuma, T. Kondo, and T. Takahashi The Europ. Phys. Jour. D, 14, 271 (2001); 

S. G. Porsev, Phys. Rev. A 56, 3535 (1997). 
[21] C. Froese Fischer, T. Brage and P. Jonsson, Computational Atomic Structure an MGHF Approach (Inst, of Phys, London, 

2000). 

[22] I. P. Grant, in Methods in Computational Chemistry, Vol. 2, edited by S. Wilson (Plenum, New York, 1988), p. 1. 

[23] F. A. Parpia, C. Froese Fischer and I. P. Grant, Comp. Phys. Comm. 94, 249 (1996); L. Visscher, O. Visser, P. J. C. Aerts, 

H. Mcrcnga, and W. C. Nieuwpoort, Comp. Phys. Commun. 81, 120 (1994); K. G. DyaU, I. P. Grant, C. T. Johnson, 
F. A. Parpia, and E. P. Plummer, Comp. Phys. Comm. 55, 425 (1989); J. P. Desclaux, D. F. Mayers, and F. O'Brien J. 
Phys. B 4, 631 ( 1971). 

[24] B. R. Judd in Atomic, Molecular and Optical Physics Handbook edited by G. W. F. Drake ( AIP, New York, 1996), p. 56 

and p. 88; T. Kagawa, Y. Honda, and S. Kiyokawa, Phys. Rev. A 44, 7092 ( 1991). 
[25] K. K. Mon, and J. B. French, Ann. Phys. (N.Y.) 95, 90 (1975); L. Benet, T. Rupp, and H. A. Weidenmiiller, Ann. Phys. 

(N.Y.) 292, 67 (2001). 

[26] A. Stuart, and J. K. Ord, Kendall's Advanced Theory of Statistics, fifth edition of Volume 1: Distribution Theory (Oxford 

University Press, New York, 1987). 
[27] Ph. Jacquod, and D. L. Shepelyansky, Phys. Rev. Lett. 79, 1837 (1997). 
[28] Ph. Jacquod, and I. Varga, Phys. Rev. Lett. 89, 134101 (2002). 
[29] I. Varga, and J. Pipek, Phys. Rev. E 68, 026202 (2003). 
[30] V. K. B. Kota, and R. Sahu, Phys. Rev. E 64, 016219 (2001). 

[31] M. Abramowtiz, I. A. Stegun (Eds.), Handbook of Mathematical functions, NBS Applied Mathematics Series, Vol. 55, ( 

U.S. Govt. Printing Office, Washington, D.C., 1964). 
[32] V. K. B. Kota, and R. Sahu, Phys. Rev. E 66, 037103 (2002); M. Horoi, V. Zelevinsky, and B. A. Brown, Phys. Rev. Lett. 

74, 5194 (1995). 

[33] V. K. B. Kota, and K. Kar, Phys. Rev. E 65, 026130 (2002). 

[34] T. Papenbrock, and H. A. Wcidcmniiller, Phys. Rev. Lett 93, 132503 ( 2004). 

[35] J. M. G. Gomez, K. Kar, V. K. B. Kota, J. Retamosa, and R. Sahu, Phys. Rev. C 64, 034305 (2001); M. Horoi, J. Kaiser, 
and V. Zelevinsky, Phys. Rev. C 67, 054309 (2003); J. M. G. Gomez, K. Kar, V. K. B. Kota, R. A. Molina and J. Retamosa, 
Phys. Rev. C 69, 057302 (2004). 

[36] V. K. B. Kota and R. Sahu, Phys. Rev. E 62, 3568 (2000). 

[37] F. S. Chang, J. B. French, and T. H. Thio, Ann. Phys. (N.Y.) 66, 137 (1971); S. S. M. Wong, Nuclear Statistical 

Spectroscopy ( Oxford University Press, New York, 1986). 
[38] G. F. Gribakin, S. Sahoo and V. Dzuba ArXive:pliysics/0401157. 
[39] A. A. Gribakina, V. V. Flambaum and G. F. Gribakin, Phys. Rev. E 52, 5667 (1995). 
[40] V. Zelevinsky, B.A. Brown, N. Frazier, and M. Horoi, Phys. Rep. 276, 85 (1996). 
[41] T. Papenbrock and H. A. Weidenmiiller, arXiv:nucl-th/0403041. 

[42] I. Lindgren and J. Morrison, Atomic Many-Body Theory 2nd ed (Springer- Verlag, Berlin, 1986). 



14 

[43] Y. Alhassid, H. A. Weidenmueller and A. Wobst arXiv:cond-mat/0406495. 



15 



State density 



CHAOS MARKERS FOR EGOE(1 +2) 
H=Ah(1)+A,{V(2)} 

(m,N,A.?t) 

Gaussian 



Fiuctuations 



Poisson 



Wigner— IDyson (GOE) 



A— r 



Strenjgtli 
functions 



5- 

function 



Breit— Wigner 



[BW domain] 



Gaussian 
[Gaussian domain] 



defines thermodynamic region 
(duaiity, different entropies coincide) 



Si 



0.08 — 



FIG. 1: Chaos markers for EG0E(l+2) 
m=6, N = 12 





oc=1.8 - 


0.1 2 
0.09 
0.06 




^ A^O. 15 - 

i I oc=2.4- : 






0.03 










O 








\ oc=3.3 

\y^=i.oi I 


0.08 

0.06 




r\ Oc=5.3 - 
/ \ Y2=0-42 : 






0.04 










0.02 

















O.OG 




jP=\ ?L=0.28 - 


0.06 


- - 






/ \ oc=9.8 - 




: A\ cx=i7' : 


0.04 




/ \ y^=0.\2.- 


0.04 


- / \ '^^^^-^ - 


0-02 






0.02 















CE— e)/CT 



FIG. 2: Strength functions for EG0E(l+2) as a function of the interaction strength A. 
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FIG. 3: Binned density of states p{E) of Nd, Pm and Sm. The dashed and dash-dot curves are the Gaussian and Edgeworth 
corrected Gaussian representations. 
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FIG. 4: (a)The energy range of CSFs for Nd arising from a non-relativistic configuration. The energies are calculated with 
respect to the lowest CSF energy within the manifold ( -9625.2394 hartree). The number of CSFs arc given within parentheses 
and the cross marks the location of the centroid. (b) The solid and dot-dash line are partial Gaussian densities p^'^{E) of 
each non-relativistic configurations and sum of all the partial densities respectively. The histogram is that of the density of 
states p"{E). 
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FIG. 5: The plots in the three rows are the averaged strength function Fk{E) for Nd, Pm and Sm respectively. As labeled in 
the plots, the plots in the four columns show Fk{E) averaged in the range —1.5 ± 0.25crfc, —1.0 ± 0.25crfc, —0.5 ± .25crfc, and 
0.0 ± .25cTfc respectively. 
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FIG. 6: The number of principal components ^2{E) = {J2^^ |Cf 1*)^^ of Nd I, Pm I and Sm I are plotted in the first row. The 
plots in the second row are binned ^2{E) and the last row is the binned density of localized states pioc(-E'). In (b) and (c) the 
continuous curves are from Eq. (24) as described in the text. 
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FIG. 7: The solid line is the number of principal components ^2{E) of energy states F around the centroid. The dot-dash 
line show the occupation number (ria)^ of selected shells which anti-correlate with the localized states. The occupation 
numbers are scaled and shifted for better comparison with $,2{E): for Nd the plot shows {n^d^^r^)^ + (nsd^^^)^ — 1 , for Pm it is 
0.2({n5d3 ,2)^-1- (nsdj ,2)^) ^^'^ ^'^^ dot-dash and dashed lines are (nes^ ,3)^ ~3 and {«6pi /j)^ + {'T-6p3,2)^ — 3 respectively. 
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FIG. 8: The \C^\'^ of neighboring localized and chaotic pair of \E) arc shown in the plots in the first column. Even though 
each pair arc neighboring states, a large difference in the structure is discernable. For the localized states, the second column 
shows the individual \k) contribution |C|^|^(na)fe to the sub-shell of highest (ria), whereas for chaotic states all the sub-shells 
have similar \C^\^ {na}k and the plots show contribution from a randomly chosen sub-shell. 
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FIG. 9: The binned Hamiltonian matrix of Nd I, Pm I and Sm I. 



